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Abstract 

This work derives the linearized equations of motion, the Lagrangian 
density, the Hamiltonian density, and the canonical angular momentum 
density for general perturbations [oc exp(im0) with m = 0, ±1, ..] of a 
geometrically thin self-gravitating, homentropic fluid disk including the 
pressure. The theory is applied to "eccentric," m = ±1 perturbations 
of a geometrically thin Keplerian disk. We find m = 1 modes at low 
frequencies relative to the Keplerian frequency. Further, it shown that 
these modes can have negative energy and negative angular momentum. 
The radial propagation of these low frequency m — 1 modes can transport 
angular momentum away from the inner region of a disk and thus increase 
the rate of mass accretion. Depending on the radial boundary conditions 
there can be discrete low-frequency, negative-energy, m = 1 modes. 
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1 Introduction 

Problems of long standing interest concern the linear modes and instabilities of 
accretion disks and the disks of spiral galaxies (for example, Kato, Fukue, & Mi- 
neshige 1998; Binney & Tremaine 1987). Interest in unstable modes in accretion 
disks has been stimulated by recent observations of quasi-periodic oscillations 
(QPOs) in accreting neutron stars and black hole candidates in binary systems 
(see review by van der Klis 2000). On a larger scale, recent observations of disk 
galaxies reveal that about 30% of the disks are "lopsided" (Baldwin, Lynden- 
Bell, & Sancisi 1980; Rix & Zaritsky 1995; Kornreich, Haynes, & Lovelace 1998). 
In some cases the lopsidedness may be dynamical in origin and an indication 
of an unstable one-armed trailing spiral wave (Lovelace, Zhang, Kornreich, & 
Haynes 1999; hereafter LZKH). Hybrid, N-body/hydrodynamic simulations of 
impulsively perturbed disk galaxies show long-lived trailing spiral features (Zelt- 
wanger et al. 2000; Kornreich et al. 2002). 
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In previous work we derived equations of motion, a Lagrangian, Hamiltonian, 
and canonical angular momentum for the eccentric perturbations [oc exp(±i(j))] 
of a thin self-gravitating disk of a galaxy consisting of a large number of rings 
(LZKH). This work utilized a Lagrangian description of the perturbation. The 
present work generalizes LZKH by deriving equations of motion, a Lagrangian 
density, Hamiltonian density, and canonical angular momentum density for gen- 
eral perturbations [oc exp(im0) with m = 0, ±1, ..] of a self-gravitating fluid 
disk including the pressure. The theory allows the calculation of the energy and 
angular momentum of the different modes and identification of modes with neg- 
ative energy and/or negative angular momentum. The negative energy modes 
can become unstable in the presence of dissipation. Negative energy modes are 
important in plasma physics (Coppi, Roscnbluth, & Sudan 1969), in Rossby 
r-mode oscillations of rapidly rotating neutron stars which are made unstable 
by gravitational radiation (Andersson 1998; Schenk et al. 2002), and for short 
wavelength modes in the shearing sheet model of modes of disks (Goodman, 
Narayan, & Goldreich 1987). We apply our theory to the "eccentric," m = 1 
perturbations of a geometrically thin Keplerian disk. We find m = 1 modes 
with frequencies much less than the Keplerian frequency. These modes can 
have negative energy and negative angular momentum. The radial propagation 
of such a mode may act to remove angular momentum from the inner region of 
a disk. Earlier, Nowak and Wagoner (1991, 1992) analyzed the modes of non- 
self-gravitating disks using a Lagrangian description of the perturbation. The 
important role of m = 1 modes for accretion disk angular momentum transport 
has been discussed by Fridman et al. (2003). Adams, Rudcn, and Shu (1989) 
studied the growth of global gravitational instabilities in star/disk systems and 
found modes which grow on a nearly dynamical time scale. Shu et al. (1990) 
presented an analytical description of the corresponding modal mechanisms. 
Tremaine (2001) analyzed slow modes using softened gravity and found these 
modes to be stable. 

In §2 we derive the linearized equations of motion and obtain the Lagrangian 
density. From this the Hamiltonian density and canonical angular momentum 
density are derived in the usual way. In §3 we introduce a representation for 
the perturbation amplitudes in terms of complex quantities. In §4 the WKBJ 
approximation is introduced for the radial dependence of the perturbations, and 
the results of the present approach are contrasted with the standard analysis. 
In §5 we consider the inner region of a low mass disk around a black hole and 
show that one of the m = 1 modes has a very low frequency (compared with 
the rotational frequency). This mode is shown to have negative energy which 
means that it may be made unstable by dissipation. In §6 we give conclusions 
of this work. 

2 Theory 

We use an inertial cylindrical (r, <f>, z) coordinate system and assume that the 
disk is geometrically thin with half-thickness h <C r. Further, we consider 
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Figure 1: Diagram illustrating equation (1), where "equil." indicates the equi- 
librium path of a fluid element and "pert." indicates the perturbed path. See 
Frieman & Rotenberg (1960). 



conditions where the accretion speed — v r <C v^. We use a Lagrangian repre- 
sentation for the perturbation as developed by Frieman and Rotenberg (1960) 
and Lynden-Bell and Ostriker (1967) and as applied to a system of rings repre- 
senting a self-gravitating disk (LZKH; Lovelace 1998). The position vector r of 
a fluid element which at t = was at r is given by 

r(t)=r (t)+€(r ,t). (1) 

That is, r (t) is the unperturbed and r(t) the perturbed orbit of a fluid element 
as sketched in Figure 1. 

Further, the perturbations are assumed to consist of small in-plane displace- 
ments of the disk matter, 

€ = £rf + &0 • (2) 

The £ r and ^ are proportional to exp(im</>) with m = 0, 1, 2, ... the azimuthal 
mode number. (The negative values of m are redundant.) 

From equation (1), we have v(r, t) = v (r ,i) + d£/dt + (v • V)£. The 
Eulerian velocity perturbation is 5v(r,i) = v(r,i) — v (r,i). Therefore, 

5v(r,t) = § + (vV)£-(€-V)v. (3) 
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The components of this equation are 



5v r = V^ r , 

Svt = V^-rfl'tr , (4) 

where 

with Q(r) the angular rotation rate of the equilibrium disk and f2' = dQ/dr. 
The equation of motion is 

^ = SF = -^\7SP+ S ^VP-VS^ , (5) 

where 6F is the perturbation force per unit mass, is the perturbation of the 
gravitational potential, £ is the disk's surface mass density, and P is the verti- 
cally integrated pressure in the disk. For simplicity we consider a homentropic 
disk where P = const S r . Under this condition, SF = — V(5P/E) — V<5$. 
We have 

d Sv d <5v . _. . _. , . 

_ = _ + (v-V)Jv+(<fvV)v. (6) 

The components of equation (6) give 



d Sv 
dt 

dSv 
dt 



= V5v r — 2Q8v ( f > , 

= (V 2 + 2Q.rO!) £ r - 2flV^ , 

= VSv^ + {K 2 r /2Q)5v r , 

= 2? 2 ^ + 2f>P£ r , (7) 



where k 2 = (l/r 3 )d(r 4: £l 2 )/dr is the radial epicyclic frequency (squared). For a 
low mass, geometrically thin Kcplerian disk, k 2 = Q 2 + 0(c 2 /r 2 ) w Q 2 . More 
generally, for ft oc l/r q , n 2 = 2(2 — q)Q 2 . For a Kcplerian disk, q = 3/2; for a 
spiral galaxy with a flat rotation curve, q = 1; and for a rigidly rotating disk, 
for example, the inner region of a spiral galaxy, q = 0. 

For SF = 0, equations (7) with £ r , ^ oc exp(im0 — iuot) can be solved 
to give four roots: two have u> — mfl = and £ r = 0, £^ ^ 0, and two have 
lu — mfi = ±n r . For the second pair of modes, TZe(£_ r ) = acos(n r t) and TZe(^) = 
— (2Q/K)asm(K r t), where a =const. and TZe indicates the real part; the fluid 
particle moves along a retrograde ellipse which is elongated in the azimuthal 
direction since n r < 2fi for most conditions (Binney & Tremaine 1987, p. 132). 
For a rigidly rotating disk, 2Q/n r = 1 so that the ellipse is a circle. 

The perturbation of the surface mass density of the disk obeys 

—jjj- + V • (£ Sv + <5S v) = , 
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where E(r) is the surface density of the equilibrium disk. Because V • (E v) =0, 
this equation implies 

5E = -V • (£ • (8) 
The perturbation of the gravitational potential is given by 

M>M) = - G /dV|^, (9) 

where the integration is over the surface area of the disk, and G is the gravita- 
tional constant. 

2.1 Perturbation Amplitudes 

We can write in general 

£r = £x cos(m0) + e y sin(m</<) , 

£;</> — —S x sin(m^) + 5 y cos(m<ft) , (10) 

where e x . y (r 7 t) and 6 XiV (r,t) are the displacement amplitudes for the 'ring' of 
the disk at radius r. For m — 1, e Xjy represents the shift of the ring's center, and 
§x,y represents in general both the shift of the ring's center and the azimuthal 
displacement of the ring matter (see LZKH). 

2.2 Equations of Motion 

From equations (5) - (7) and (10), we obtain the disk equations of motion, 

e x + 2mfle y — m 2 Q 2 e x — 2Q(5 y — mQ,S x ) = 5F^ , 

e v - 2mCle x - m 2 n 2 e y + 2fl(S x + mQ8 y ) = SF^, 

S x + 2mflS y - m 2 n 2 S x - 2fl(e y - mfle x ) = -5F^, 

5 y -2m,n5 x -m 2 ri 2 5 y + 2n(e x + mne y ) = 5F^, 



(11) 



where e x = de x (r,t)/dt, etc. Here, 

n 2 = fl 2 - 2VtrQ! /m 2 , 

and 



6F^ S = l^j- [cos(m<£), sin(m^)] SF a 



with the S, C superscripts denoting the sine and cosine components, and a = 
r, (j). 

The force per unit mass 



SP 

SF a = -V— - VS<S> = V 

2-i 



c 2 s V • (££) 



V(5$ , (12) 
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with a = r, <j>, and c 2 — dP/dYj is the effective sound speed (squared). The 
pressure force on the right-hand side of equation (12) gives 

= D[c 2 s D^e x )/T]-D{clm5 x /r) , 

SF S / = D[clD^e y )m-D(c 2 s m5 y /r) , 

-SF S / = [mc 2 /{r^)]D^e x )-c 2 s m 2 5 x /r 2 , 

5Ff = l m c 2 /(rX)}D^e y )-clm 2 5 y /r 2 , (13) 



where £>(..) = d(..)/dr and 


= (l/r)d(r ..)/dr. 


The gravitational force gives 






6F? 9 = 


-lfr G J 


f d 2 r' K m R' x 


5F r S9 = 


dr G J 


f d 2 r' JC m R' y 


-6F*° = 


-7°) 


1 d 2 r' K m R' x 


5F°° = 


~r G .1 


f d 2 r' K m R' y 



where d 2 r' = 2irr'dr', R'. = Rj{r',t), 



(2m- 1)!! yp— (7)2j 



(14) 



3 ~ "J* 

Rj{r) = D*(£e,) - mUSj/r (15) 

for j = x, y, and 

(r, r ) = — * . (16) 

2tt 7 ^r 2 + (r') 2 - 2rr' cos(l') 

We can express this as 

( _ 2)m 



where 2 = [r 2 + (r') 2 ]/|r 2 - (r') 2 \, (-1)!! = 1, and P™ is the usual associated 
Legendre function. 

2.3 Lagrangian 

The Lagrangian which gives the above equations of motion is 

L = J d 2 r £ = 2ir J rdr C , (17) 
where £ = Cq + L\ + Li + £3 is the Lagrangian density with 



e 2 + <5 2 + m 2 2 e 2 + m 2 2 <5 2 - 4mf> 2 e • 5 



G 



Ci = -SO me x e + mS xS-Sxe-exS ■ z 



c 2 = -i|[v.(so] 2 

1„ , m<5 x 2 



£3 = f E (Q - ^) ■/ rf 2 r* m (r./)E' ( Q < - ^) . (18) 

Here, e = e x x + e y y, S = S x x + S v y, Qj(r) = D^He^/T, (j = x,y), Q = 
QccX + QyY, and £' = £(/). 

The Lagrangian density is C = £[e, S, e, 5, Q]. In terms of functional deriva- 
tives, the equations of motion are 

d SL SL d ( 1 SL 

■ 



dt 8ij Sej dr \S SQ 

= *± (19) 

(see for example Goldstein 1950, eh. 11). For the contribution from £3, 
5L 3 /SQ = GE / d 2 r'JC m E'(Q' -mS' /r') and SL 3 /SS = -G£(m/r) / dV/C m S'(Q'- 
md'/r'). 

2.4 Hamiltonian 

The Hamiltonian density is 

«-*■%+*■%-'■ m 

The time derivative of this equation gives 

1 d . r , „ . <9£ 

0F + r* (r ^ ) = o ' ^ e -aq< (21) 

where J-^r is the radial energy flux-density. The explicit Hamiltonian density 
is 

1^9 2 



+ 2^{ Q - r 
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in terms of the Lagrangian variables. Assuming vanishing energy fluxes through 
the inner and outer radii of the disk, we have 

H = J d 2 rU , with ^ = , (23) 

where H = E is the energy of the perturbation. 

2.5 Canonical Angular Momentum 

We can make a canonical transformation from rectangular to polar coordi- 
nates, (e x ,e y ) — > (e, ip), (6 x ,S y ) — > (S, ip) to obtain the Lagrangian as C = 
C(e, 5, p, ip, e, 5, ip, ip, D*(Se)/S, dp/dr). Note for example that e 2 — > e 2 + e 2 ^ 2 . 
We then find that C is invariant under the simultaneous changes ip — > <y9 + 0, 
ip ^ ip + 9. Therefore, it is useful to make a further canonical transformation, 
[y, V'] ~ > i a = (f + "0)/2,/3 — y — ^>], which gives dC/da = 0. Defining the 
canonical angular momentum density as 

we obtain 

dVd, 1 9 , ^ , „ ^ dC 
ot r or o(oa/dr) 

where T-p r is the radial flux-density of angular momentum. Assuming vanishing 
angular momentum fluxes through the inner and outer radii of the disk, we have 

d 2 r V<p , with -± = , (25) 

where is the canonical angular momentum of the perturbation. From equa- 
tion (20) we find 



e 2 {p - mfi) + S 2 (ip - mfi) + 2fie<5cos(< / s - -0) 



(26) 



where the last term can also be written as 2Y,Qe • 8. 



3 Complex Displacement Amplitudes 

It is useful for later work to introduce the complex displacement amplitudes, 

£ = e x - ie y = e(r, t) cxp [ - ip(r, t)] , 

A = S x -i6 y = S(r, t) exp [ - iip(r, t)] , (27) 

(LZKH). Here, e and 5 are clearly non-negative, and the angles <p and ip are the 
line-of-nodes angles for the £ and A components of the perturbation. In terms 
of these complex amplitudes, equations (11) can be combined to give 

£ + 2imfl£ - m 2 n 2 £ - 2ittA + 2m0 2 A = SF £ , 



8 



A + 2im£lA - m 2 n 2 A - 2i£l£ + 2mtt 2 £ = 8F A , (28) 
where SF £ = 8F° - i5F s r and 6F A = -SF^ - iSF°. 

4 WKBJ Approximation 

For perturbations with short radial wavelengths A r = 27r/|fc r |, the WKBJ ap- 
proximation is useful. That is, for h < A r /(2-7r) < r, we can take £, A <~ 
exp[i J r dr'k r (r')] in equation (18). The radial derivatives can be approximated 
as for example Z?*(X£) « ik r 'E£. The azimuthal wavenumber k$ = m/r is as- 
sumed to be less than or of the order of \k r \ which corresponds to < m < r/h. 
Note that the radial wavenumber k r may be > or < 0. For m > 1, k r > 
corresponds to a trailing spiral wave and k r < to a leading spiral wave. 
With the WKBJ approximation, the pressure force becomes 

8F% = -k 2 .c 2 £ - ikrk^A , 

5F P A = -k\c 2 s A + ik r k^c\£ . (29) 

The gravitational force is 

SF§ = 2irGY,(k 2 £ + ik r k 4> A)/k , 

5F A = 27rGS(fc^A - iKk^jk , (30) 

where k = (k 2 + k 2 ,) 1 ! 2 . 

In terms of £ and A, the Hamiltonian density in the WKBJ approximation 
becomes 



H = ^(\£\ 2 + \a\ 2 -(m 2 n 2 -2nrn')\£< 2 



2 

- m 2 n 2 \A\ 2 +4m{l 2 Ke(£A*) 



V 7 ,2 

+ — \ik r £-k^A\ 2 , (31) 

where u 2 = c 2 — 2-nGY>/k, lZe(..) indicates the real part, and (..)* the complex 
conjugate. The radial energy flux-density is 

J-Er = — £ u 2 e • (Q — k^S) 

= -S u 2 Ke[£* (ik r £ ~ k.pA)] , (32) 

from equation (23). 

Similarly, the canonical angular momentum density can be written as 

= X^lm[£(£* - imn£*)] + lm[A(A* - imttA*)] 

+ 2mZe(£A*) \ , (33) 
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where Im(..) indicates the imaginary part. The radial flux-density of canonical 
angular momentum is 



T Vr = £ u 2 Jm[£* (ik r £ - fc A) ] , 



from equation (26) . 



4.1 Mode Frequencies 

Using equations (28) and taking £, A ~ exp(— iut), we obtain 





A 12 ' 




' £ ' 


_ A 21 


A 22 




A 



where 



An 
A12 
A21 
A22 



= 



-Cj 2 + 2fM' + (k r c s ) 2 - 2nGZk 2 /k 
—251a) + ik r k<pc 2 — 2jTiGT,k r k ( f ) / k 
— 2Q(D — ikrk^c 2 . + 2 / niGT,k r k ( f ) /k , 
-Co 2 + {k^Csf - 2irGZk 2 /k , 



(34) 



(35) 



(36) 



where (D(r) = a; — mfi(r) is the wave frequency in the reference frame comoving 
with the disk at r. 

The dispersion relation implied by equation (35), det(A) = 0, gives 

wi 2.3,4 = mVt ± A 



K 2 + k 2 U 2 ± B 



^/{k 2 + k 2 u 2 ) 2 - 8k 2 u 2 nrn< 



(37) 



where the (A, B) subscripts on the ± signs indicate that they are independent, 
and k 2 = k 2 + k\. We adopt the convention that uj a (A,B) = u>\(+, +), 
u)2{+, —), W3(— , +), and a>4(— , — ). The frequencies are invariant under k r — > 
— k r owing to the time reversibility of the basic equations; trailing spiral waves 
(k r > 0) are equivalent to leading waves (k r < 0). 

The mode frequencies found here differ in some respects from the standard 
analysis owing to our different application of the WKBJ approximation. In the 
standard approach, one takes, for example, V • 5v w ik • 5~v (Binncy & Tremaine 
1987, ch. 6; Kato et al. 1998, ch. 13), whereas we take V • £ « ik ■ £. Note that 
the standard approach does not lead to equations of motion derivable from a 
Lagrangian. 

For axisymmetric k^ = m/r = perturbations, equation (37) gives two zero 
frequency modes and the two modes with 



Ld 



= ±( n 2 r + k 2 c 2 -2irGZ\k r 



1/2 



(38) 
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which is the same as the classical result of Toomre (1964) and Safronov (1960). 
As a function of k r , to 2 has a minimum at krs = 7rG£/c 2 at which point 
w 2 (£; r o) = — (7rGE/c s ) 2 . If this minimum value is less than zero, which 
occurs for S > Y*ts = n r c s / (ttG) , then one of the oj values is positive imaginary 
corresponding to a linear Jeans-type instability (Safronov 1960; Toomre 1964). 
With Q = Ets/S, there is instability for Q < 1 (Toomre 1964). 

For tightly wrapped spiral waves (k 2 ^> fc 2 ,), the two modes with ±b = +1 

1 /2 

simplify to w = mfl ± [k 2 + k 2 c 2 — 27rGS|fc r |] which is the same as the 
dispersion relation of Lin and Shu (1966) and Kalnajs (1965) with the "reduction 
factor" T =\. 

Consider now equation (37) for general (fc r , k^) with |k| increasing but with 
tan(a) = \k^/k r \ =const. The argument of the inner square root in equation 
(37), (K 2 +k 2 u 2 ) 2 — 8k 2 u 2 Qrfl' , has a minimum as a function of |k| at |k| = k TS , 
the same as for k$ = 0. If the minimum of this argument is less than zero, one 
of the mode frequencies has a positive imaginary part (as well as a real part for 
oj — mCl). This condition for instability can be expressed as 



Xcrit 2qS 2 f qS 2 g 2 S 4 

Sts I 2-q [2-q + (2 - qf 



1/2' 



1/2 

» , (39) 



where S = sin(a) = |fc^|/|k| and q = — dlnfl/dlnr. For S = 1 this gives 
S c „ t /S T5 > 0.268 (or Q < 3.73) for q = 3/2 and S crit /E TS > 0.414 (or 
Q < 2.42) for q = 1. Thus, the non-axisymmetric modes are unstable at larger 
Q than the axisymmetric modes. 

For |k| > 2krs, we have k 2 u 2 > so that the argument of the inner square 
root in equation (37) is greater than zero. As a result, the outer square root in 
equation (37) for ±b = — 1 is purely imaginary if q > 0. Thus the mode with 
±,4 = +1 is unstable for q > 0. 

Figure 2 shows the four mode frequencies as a function of fc^ for a case with 

E cr ii < S < StS- 

4.2 Shearing of Unstable Modes with k^r 3> 1 

For an unstable wave with k^r = m » 1, the geometrical optics relations for a 
wave- packet give 

dk r duo , ^ d(rkd) du „ , 
" -*^< dt = ' (40) 

so that k r w qfc^rit+const., and rfc^ =const. Thus, a wave-packet initially with 
long radial wavelength, with \k r \h <C 1, is rapidly sheared by the differential 
rotation of the disk if q > 0. In a time of order t s } lear — (qk^hfl)^ 1 the wave 
evolves into a very short radial wavelength, \k r \h > 1, trailing spiral wave. 
This is discussed in detail by Goldreich and Lynden-Bell (1965). The maxi- 
mum amplification a wave can have (without reflections) is A = cxp(J dtuii) <~ 
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Figure 2: Real and imaginary parts of the four mode frequencies (labeled 1, ..,4) 
normalized by £1 for » k%, c s /(rQ) — 0.05, and S = 0.4Sts, where Sts is 
the surface density at which the Toomre-Safronov instability sets in for axisym- 
metric perturbations. We have taken q = 3/2 which corresponds to a Keplerian 
disk. For these conditions rkrs = 8. 
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exp[m.ax(u>i)t s hear]y where u>i = Xm{uj) is the growth rate. For the above- 
mentioned unstable modes for ^ cr u < S < Sts, which have k^r ^> 1, we find 
A ~ cxp(l), which is probably not significant in the absence of reflections or 
trapping of the waves in the vicinity of a bump in the disk as a function of r 
(Lovelace et al. 1999a; Li et al. 2000). 



5 Negative Energy Modes 

We consider the case of a low mass disk around a Schwarzschild black hole where 
the inner radius of the disk is at Ti — 3rs with rs = 2GM/c 2 the Schwarzschild 
radius. We adopt the pseudo-Newtonian potential $ = —GM/(r—rs) (Paczyiiski 

6 Wiita 1980; Kato et al. 1998, ch. 2). For this potential, 



n 



r<i\/r \ V2 / q \ 1/2 

GM \ n r I r — Srg x 



r(r — rs) 2 ) £1 \ r — rs 

din SI 3 r — rs/3 
^ d In r 2 r — rs 



(41) 



Figure 3 shows the radial variations of these quantities near the inner edge 
of the disk. In the following the mode frequency u> is measured in units of 
Sis = \/GM/r 3 s , while ft is in units of s _1 . Radial distance R and l/k r are 
measured in units of rs, while £ and A are in units of length. Note that 
n s w 7.16 x 10 4 (M Q /M)s- 1 and r s « 2.96 x 10 5 (M/M Q )cm. 



5.1 One- Armed Spiral Waves (m = 1) 

Figure 4 shows the radial wavenumber dependence of the low-frequency mode 
loi and the corresponding mode energy (per unit area) dE/dA (from equation 
(31)) for a sample case for k r > which corresponds to a trailing spiral wave. 
(The 1— subscript on lo is dropped in the following.) The radial phase velocity 
is outward (u/k r > 0) and the azimuthal phase velocity is positive (ujr > 0). 
[Note that both quantities plotted in Figure 4 are even functions of k r , so that 
the plot is the same for k r < which corresponds to a leading spiral wave.] 
For k r less than a certain value k r o, the wave frequency is positive while the 
wave energy is negative. Further, the energy flux-density of the perturbation 
Tet (from equation (32)) is positive (outward) for k r < k r0 and negative for 

From Figure 4 the group velocity v gr = dw/dk r is negative for k r > 0. Thus, 

the sign of TetI (dE/dA) < is the same as that of v gr for k r >k r0 as it should 
be. Note that the radial phase and group velocities have opposite signs. 

The canonical angular momentum (per unit area) V$ (from equation 33) is 
negative and the flux of angular momentum T-pr (from equation 34) is positive 

for the low frequency mode for k r >k r0 . Thus T-p r jV^, < has the same sign as 
v gr . 
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0.1 




-0.1 r 1 

1 2,3 4 5 

k r r 

Figure 4: Dependence of the low-frequency mode u> on radial wavenumber from 
equation (37) and the corresponding wave energy (per unit area) dE/dA from 
equation (31). For k r < note that both curves are even functions of k r . Here, 
us is normalized by Cls, and dE/dA by Ef2 2 |£| 2 /2. For this plot R = 10 and 
c s /(rVL) = 0.1. 
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Tabic 1: Physical Quantities for the Low- Frequency Mode 





uj dE/dA 


J^Er v gr = duj/dk r Vtf, 




< k r < k r0 


+ 


+ 


+ 




+ 




+ 



The signs of the different quantities are summarized in Table 1 for the trailing 
waves (k r > 0). For the leading waves, k r — > — k ri the signs of dE/dA and 
do not change, but the signs of the fluxes Tet and Fvr reverse as does the sign 
of the group velocity v gr . 

In contrast, for the high-frequency mode u>2, the wave energy is positive 
independent of k r . Also, the angular momentum is positive. 

For the low-frequency mode the approximate dependence is 

— ('"IS)- 

Figure 5 shows the R— dependences of cj a and (k r r) for c s /(rQ) = 0.1. For 
this case uoq « C^/R 11 with C w = 4.71, a = 3.16 and (k r r) w Ck/R b with 
Cfe = 21.1, b = 0.702. The frequency is of course independent of R so that 

^ = lS,fl-^ (43) 



R 2 + 2b \ a 

where the allowed region k 2 > extends from the inner radius of the disk 
R = Ri = 3 out to the turning point R tU rn = (C LU /uj) 1 ^ a > 3. 
Because dE/dA w -w£ri 2 |£| 2 /2, we have 

/•OO 

E^-irr 2 s uj RdRW 2 \£\ 2 , (44) 
J ft, 

which is negative in that u) > for the 'bound states'. Also, we find V^, ~ 
— Sri|£"| 2 so that the total canonical angular momentum of the disk perturba- 
tions is 

/•oo 

L z » -tit! / RdR £0|£| 2 . (45) 
Jft, 

We discuss the mode energy and angular momentum further in §5.3. 

In the limit u — > 0, the trailing spiral wave for example has the form R = 
R l /(l + 63V/C fe ) 1/b for (j> < 0. 

5.2 Standing Waves 

A standing wave solution for say £ can be written as 

£ = £+ + £-, 

= A|fc r |- 1 / 2 exp(^) + B\k r \^/ 2 cxp(-i-f) , (46) 
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3 4 5 6 7 8 910 20 

R 



Figure 5: Radial dependences of the two 'constants' in equation (42) for the 
case where c s /(rO) = 0.1. The straight lines correspond to least squares fits 
through the shown points which give (fc r r) = 21.1/i? 702 and uj = 4.71/i? 316 . 
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where the £ + term corresponds to a trailing spiral wave and the £_ term a 
leading wave. Also, 



is the phase with the u> dependence made explicit. An expression for A in terms 
of £ can be derived from equations (35) and (36). 

The specific nature of the modes clearly requires a definite boundary con- 
dition at the inner edge of the disk. The present model, which neglects v r , 
does not allow a full treatment of the inner radius of the disk where the radial 
inflow may be transonic (see for example Kato et al. 1998, ch. 9). Instead, 
we first consider an inner disk boundary condition appropriate for the Shakura 
and Sunyaev (1973) model where for r — > r i; £ — > oo; namely, £{ji = 3r^) = 0, 
which corresponds to perfect reflection. For this boundary condition, and the 
usual WKBJ matching at the outer turning point Rtum, we obtain the Bohr- 
Sommcrfcld quantization condition 



which determines the allowed ui values. The 3/4 in equation (48) (rather than 
1/2) is due to the boundary condition at Ri. For the mentioned dependences of 
wo and (k r r)o, and c s /(rO) = const, there are a finite number of modes n with 
the lowest frequency corresponding to the largest n value. 

Instead of using the approximate equation (48), we solve the wave equation 
for the perturbation amplitude, 



This is done using a shooting method: We guess u, integrate outward from 
R = Ri = 3, and require £{R — > oo) — > 0. This requirement determines the 
allowed frequency values. Figure 6 shows the behavior of the wave functions for 
a sample case with c s /(rQ) = 0.1. The corresponding frequencies in units of Qs 
are uj n w 3.75 x 10~ 2 , 7.74 x 10~ 3 , and 8.69 x 10~ 4 for n = 0, 1, 2, respectively. 
For a M = 1OM black hole, these frequencies in Hz are: 268, 55.4, and 6.22. 
These frequencies agree approximately with those obtained from equation (48). 
This equation implies that there are only three modes (or bound states) for 
c s /(rQ) = 0.1. 

Figure 7 shows the variation of the mode frequencies with the disk thermal 
speed c s /(rO). 

With regard to the energy of the standing wave modes, notice that we now 
have in place of equations (28) separate equations for £+, £-, A+, and A_. 
Thus the Lagrangian density becomes C = £+(£+,..) + £_(£_,..), and the 
Hamiltonian H = H++H-. The energy of the trailing and leading components 
are equal and negative. Thus the total energy E = E + + E- is negative. From 




(47) 




n = 0, 1,2, .. 



(48) 




(49) 
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Figure 6: Wave function profiles £ n (R) for the three low- frequency standing 
wave modes for the case c s /(rtt) = 0.1 obtained by solving equation (49) using 
a shooting method. The normalization of the different wavefunctions results 
from the boundary condition (d£ n /dR)i = 1. 
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1 




0.25 



c s / (rfl) 



Figure 7: Dependence of mode frequencies uj„ in units of fls on the disk temper- 
ature as measured by c s /(rCl). The curves are smooth fits through the shown 
points. 
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270 



Figure 8: Time dependence of the linc-of- nodes angle ip(R,T), which is such 
that £ = \£{R, T)\ exp[— iip(R, T)], obtained by solving the Schrodingcr equation 
(51) for an initial wavefunction £(R,T = 0) = const(i? — 3) cxp(— i?/4) which 
is not an eigenfunction. The boundary conditions were £{Ri) = = £{R ma x) 
and a nonuniform mesh Rj = RiQi , j = 1, ..,300 with Q — 1.009 so that 
Rmax ~ 44.1. The time-dependent evolution represents a superposition of the 
three eigenfunctions shown in Figure 6 with the average of the energy < E > = 
—8.09 x 10~ 3 w. The plot shows the reflection of a trailing spiral wave (T < 125) 
from the inner disk radius which gives rise to a leading spiral wave (T > 125) 
for 3 < R < 5.6. 
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equation (32) the total energy flux is 



T Er = w ^Im(^+£! r 



dR - dR 



dR, 

= ^|fc r |(|£ + | 2 -|£_| 2 ) . (50) 

For the considered boundary condition at the inner disk radius R = 3, Tet = 0. 
At large distances, Tet = because £(R —> oo) = 0. Thus from equation (21) 
we have dE/dt = 2n(rJ r Er)i — 27r(r^ 7 E r ) 00 = 0, where the i— subscript indi- 
cates evaluation at R = 3. Similar conclusions apply to the canonical angular 
momentum where T-p r = £ c slM(|£+| 2 — l^-| 2 )- 

5.3 Time-Dependent Solutions 

Equation (49) can be written in the form of a Schrodinger equation 

d 2 



w d£ 



dR 2 



U(R) 



£ , (51) 



where T = Sl s t, U(R) = -C 2 /R 2+2b andw = (C 2 / 'C u ) / 'R 2 + 2l - a . The R- dependence 
of w is rather weak (~ i? -0 24 for c s j (rO) = 0.1) and for this reason we make the 
simplification that w = w(R = R4). The sign of the time derivative is opposite 
that of Schrodinger equation because the energy is negative for oj > 0. 

For the mentioned bound states, E n = — wuj n , where the tilde indicates 
a different normalization of the energy compared with equation (44). For the 
unbound states or continuum modes we have LOk < and Ek = — w Uk > 0. For 
the continuum trailing waves (k r > 0) both the radial phase and group velocities 
arc negative, whereas for the leading waves (k r < 0) both are positive. 

Figure 8 shows polar plot (see Briggs 1990; LZKH) of the time-dependence 
of the mode for an assumed initial perturbation £(R,T = 0) = const(_R — 
Ri) exp(— -R/4) which is neither leading nor trailing. The mode evolves to a 
predominantly trailing spiral wave. Trailing spiral one-armed spiral waves were 
found under different conditions in our earlier study of galactic disks (LZKH). 



6 Conclusions 

The perturbations of thin disks have been studied in many works applied to 
understanding the waves or modes in galactic disks and in accretion disks. This 
work has derived the linearized equations of motion, the Lagrangian density, the 
Hamiltonian density, and the canonical angular momentum density for general 
perturbations [oc exp(zm</>) with m = 0, ±1, ..] of a geometrically thin self- 
gravitating, homentropic fluid disk including the pressure. The theory was ap- 
plied to the "eccentric," m = ±1 perturbations of a geometrically thin Keplerian 
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disk. The m — 1 modes may be at low frequencies compared to the Keplerian 
frequency. Furthermore, the modes can have negative energy and negative an- 
gular momentum. Their propagation can remove angular momentum from the 
inner region of a disk and thus act to enhance the accretion of matter. De- 
pendent on the radial boundary conditions, there are discrete, low-frequency, 
negative energy m = 1 modes. Time-dependent perturbations tend to give rise 
to predominantly trailing spiral waves. A future direction for extending this 
paper include the treatment of the non-linear evolution following the approach 
of Shu, Yuan, and Lissauer (1985), Heemskcrk, Papaloizou, and Savonije (1992), 
and Lee and Goodman 1999). 

We thank M.P. Haynes, D.A. Kornreich, N.F. Comins, and M.M. Romanova 
for stimulating discussions on eccentric perturbations of disks. This work was 
supported in part by NASA grants NAGS-9047, NAGS 9735, by NSF grant 
AST-9986936, and by CRDF grant KP-2555-AL-03. 
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